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Abstract: Supersymmetric extensions of the Standard Model with highly mixed 
squark flavours beyond minimal flavour violation provide interesting scenarios of new 
physics, which have so far received limited attention. We propose a calculable re¬ 
alization of such scenarios in models of gauge mediation augmented with an extra 
interaction between the messengers and the up type squark. We compute the super- 
symmetric spectrum and analyze the flavour physics constraints on such models. In 
a simplihed model approach, we show that scenarios with maximal squark flavour 
mixing result in interesting phenomenological signatures at the LHC. We show that 
the model can be probed up to masses of rriu < 950 GeV in the single-top event 
topology at LHC14 with as little as 300 fb“^. The most distinctive signature of 
highly mixed scenarios, the same sign positive charge di-top, can be also probed to 
comparable squark masses at high luminosity LHC14. 
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1 Introduction 

Recent results from the Large Hadron Collider (LHC) impose stringent limits on 
the scale of supersymmetry and considerably constrain the minimal supersymmetric 
extension of the Standard Model (MSSM). Yet, the spectrum of non-standard su¬ 
persymmetric (SUSY) models has not been studied in full detail, while the current 
data constraints on non-standard SUSY are weaker [1-6]. At the dawn of the LHC 
Run H at = 14 TeV, and with the prospects for the future high luminosity LHC, 
it is important to revisit constraints imposed in studies of the MSSM parameter 
space, and to explore if some viable realizations have been overlooked and/or could 
be accessed through new types of collider signatures. 

As an example of scenarios which were previously less studied, one could consider 
relaxing the minimal flavour violation (MFV) assumption on SUSY models. One 
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interesting option for non-standard SUSY beyond MFV involves considerations of 
mixing between sqnark flavonrs. The immediate benefit of scenarios where sqnark 
mixing is allowed is that branching ratios of sqnark decays into hnal states predicted 
by standard SUSY models are lowered, hence weakening the experimental bonnds 
on the stop mass. Cnrrent experimental resnlts, particnlarly from flavonr physics, 
give strong constraints on mixing between hrst and second generation sqnarks, while 
mixing between third and hrst/second generations remains relatively nnconstrained. 
Several scenarios of SUSY models with sqnark mixing have recently been discussed 
in Refs. [6-10]. 

In this paper we investigate the possibility that the low energy part of the super- 
symmetric spectrum is characterized by a single (right handed) light sqnark, highly 
mixed in havour between stop and up-squark or stop and scharm-squark. We pro¬ 
pose an explicit realization of such scenarios in a simple extension of gauge mediated 
models, augmented with an extra interaction between the messengers and the right 
up type quark superheld. We study the constraint imposed on such models by havour 
observables and by the requirement of a viable Higgs mass. A similar extension of 
the gauge mediation scheme has been recently studied in several diherent contexts 
[11-27], for instance with the purposes of addressing the problem, of raising 

the Higgs mass and/or investigating their possibly unusual havor patterns. 

We study in detail the reach of LHC14 for the case where the light sqnark is max¬ 
imally mixed between top- and up-type sqnark. We focus on the collider signatures 
which are peculiar to a scenario with maximal stop-sup mixing, i.e. single top and 
same sign di-top production. These signatures have already been investigated in sev¬ 
eral new physics scenarios. The single top signatures have been studied for instance 
in [28-33]. Same sign tops have been studied in new physics models attempting to 
explain the forward backward asymmetry [34-40], or in composite models [41, 42]. In 
the context of supersymmetry, same-sign top hnal states can arise either in R-parity 
violating models [43-45] or in R-parity preserving theories through a gluino decay 
chain [46, 47]. However, in the latter case we expect a roughly equivalent number of 
same sign positive top pairs and of same sign negative top pairs. Instead, in scenar¬ 
ios with large sup-stop mixing, the same-sign top contribution at the LHC will be 
dominantly of positive charge, since it is obtained through an initial state of two up 
quarks. In the context of R-parity preserving MSSM, the presence of the single top 
signature in association with the same sign positive top signature would hence be a 
robust hint of large sqnark mixing. 

We dehne a simplihed model for LHC searches that consists only of the gluino 
(that we hx representatively at 2 TeV), the maximally mixed sup-stop right-type 
sqnark and the neutralino (mostly Bino). Considering benchmark points not yet ex¬ 
cluded by LHC8, we show that the single top signal can be discovered at LHC14 with 
300 fb“^ up to rriu-^ < 950 GeV in the aforementioned channels. The more distinc¬ 
tive signature of sizeable sup-stop mixing, (i.e. same sign tops), leading eventually 
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to same sign positive leptons plus missing energy, can be probed at High Luminos¬ 
ity LHC 14. Our collider analysis represents a concrete proposal to test the mixing 
property of the light squark in the next LHC runs. 

The paper is organized as follows: In section 2 we present the model setup and 
we explore the relevant parameter space, identifying the phenomenologically viable 
regions and possible benchmark points. On the basis of this analysis, in section 3 we 
define a simplified model for collider studies, with a light and maximally mixed sup- 
stop squark (MMUT). In section 4 we discuss its main production modes at collider 
and the constraint on the model from LHC8. In section 5 we study in detail the 
LHC14 signals of the simplihed model MMUT, i.e. processes with a single top quark 
and missing energy Ipj, or the distinctive same sign positive top with Jprp signature. 
We show the reach of LHC14 for these hnal states on some representative benchmark 
points of the MMUT model. 

2 The model formulation 

2.1 Low Energy Constraints 

In the minimal supersymmetric extension of the Standard Model (MSSM), the super- 
symmetry breaking parameters (the soft terms) can be generic sources of flavour and 
CP violation. Hence, low energy observables put strong constraints on the structure 
of the soft terms and more generally on the mechanisms of supersymmetry breaking 
and its mediation to the MSSM. In this context, gauge mediation (see e.g. [48]) is one 
of the most interesting possibilities, since it provides a predictive and computable 
framework while accommodating for the constraints from flavour physics. In this 
paper we consider an extension of gauge mediation involving extra interactions that 
break the flavour degeneracy in a controllable way, by modifying the structure of the 
soft terms. Our purpose is to provide a predictive framework for supersymmetric 
models with large mixing in the right squark up sector, and with one light squark 
eigenstate. In particular we will focus on RR mixing between the first and the third 
generation. 

The flavour bounds in the MSSM can be analysed in a model independent way 
by constraining the structure of the supersymmetry breaking parameters. Here we 
briefly review the constraints on the up-type squark mass matrix that are relevant 
in our analysis, while for a more comprehensive review we refer the reader to [49]. 

The up-type squark mass matrix in the superCKM basis can be written as 




( 2 . 1 ) 


where the dimensionless 6^ parameterise the deviations from flavour alignment. We 
are interested in flavour mixing showing up in the right-right part of the up-type 
squark mass matrix, since they are less constrained by flavour physics [2]. 
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The matrix is a 3 x 3 matrix which is determined by the soft mass of the 

up-type squark The most relevant constraint on the mass matrix is obtained 

through the measurement of the Do — Dq mixing, which bounds the absolute value 
of the to be smaller than 0.05 [49, 50]. Apart from this bound, low energy 

observables do not independently constrain 13 or (5^^)23- 

Another possibly relevant constraint to keep into account for the 1-3 mixing is 
given by the neutron electric dipole moment (EDM). Even in the case of no CP phases 
in the supersymmetry breaking parameters, a small deviation of the soft terms from 
being proportional to the identity matrix can induce some extra CP violating effect 
once we rotate into the CKM basis. The strict experimental bounds on the neutron 
EDM can lead to stringent constraints on the allowed mass spectra [51]. However, 
the neutron EDM processes also involve L-R mixing, hence they can be relevant only 
in the combined presence of large 1-3 RR mixing and (diagonal or off diagonal) large 
LR mixing. 

In the following we provide a computable model which induces large off diagonal 
contribution only to the RR up-type squark mass matrix (m?)ij in the framework of 
extension of gauge mediation, and which is compatible with flavour constraints. 

2.2 Extended Gauge Mediation 

We proceed to dehne a class of models with large squark mixing in the context of 
gauge mediated supersymmetry breaking. Scenarios with large squark mixing can 
be accommodated by extensions of gauge mediated models augmented with extra 
superpotential interactions among the messengers and some superheld of the MSSM. 
The typical realization of gauge mediation (for a review see [48]) includes messenger 
helds charged under the SM gauge groups. It is then natural to explore the possibility 
that the messengers could couple directly to SM matter helds. Extending gauge 
mediation with extra superpotential interactions in turn induces extra contributions 
to the soft terms. 

This class of models has recently been studied in several papers [11-14, 16-27]. 
An interesting aspect of the above mentioned superpotential deformations is that 
they are not automatically diagonal in havour space, and hence they can represent 
controllable sources of havour violation, which are normally absent in standard gauge 
mediation.^ In addition, the set of possible extra interactions involving the messen¬ 
gers (in complete representation of unihed gauge group) and the SM matter helds 
have been classihed in [ 11 ], where also the complete formulas for the induced soft 
masses have been computed. 

Here we focus on a specihc superpotential interaction involving messenger helds 
and the right handed type squarks, since we aim for the possibility of inducing large 

^Another interesting aspect of these deformations, which we do not exploit here, is that they 
could also induce large A-terms. 
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mixing in the right up squark mass matrix. For this reason our model differs from 
the one considered in [11], since it has a generic flavor pattern. Moreover, contrary 
to most models considered in the literature, it also respects a discrete R symmetry. 
This motivates the absence of other MSSM-messenger couplings and avoids the prob¬ 
lematic issue of generating large off-diagonal A-terms which are generically strongly 
bounded by low energy observables. Large off-diagonal A-terms could also lead to 
other problematic issues in the evolution of the SUSY spectrum, as we will discuss 
in the following. 

The model we consider consists of a pair of messengers in the 5-1-5 that we denote 
{4>i,(p2,4’i,4>2)- We assume that the component of the 5 messengers with the same 
quantum numbers of right handed down quark interact with the up type quarks via 
the superpotential couplings 

3 

(5hF = A Cif/i0i(5^£))02(5,_D)) (2.2) 

i=l 

where the index i runs here on the flavour index. Hence the interaction is not diagonal 
in flavor, and can induce non-trivial flavor mixing, depending on the values of the 
vector c = { 01 , 02 , 03 }. We take the normalization \/cf~+~cf+~^ = 1, while A sets 
the overall size of the deformation. We take all couplings to be real, in order to not 
introduce sources of CP violation. The superpotential (2.2) adds to the usual MSSM 
superpotential which includes the Yukawa couplings and the /i term. 

The supersymmetry breaking superpotential for the two pairs of 5 -|- 5 messengers 
is 

hFfi.pres = + 02<^2) + U(0102 + 020l), (2-3) 

with the SUSY breaking spurion Y = 9‘^F. The messenger scale M sets the energy 
scale where the soft masses are generated. Note that there is a residual Z 4 R- 
symmetry under which the spurion Y has charge 2, with the other charges reported 
in Table 1. This implies that gaugino masses and also A-terms cannot be generated 
by this sector and by the deformation A. 



Y 

01 

01 

02 

02 

U,D,E 

Q-) L 

Z 4 R-symmetry 

2 

2 

0 

0 

2 

0 

1 


Table 1: Z4 i?-symmetry charge assignment. 


The discrete R-symmetry, together with a messenger Z 2 symmetry under which 
and are odd, implies that the deformation ( 2 . 2 ) is the only one compatible with 
these discrete simmetries and with gauge invariance. Indeed the Z 2 symmetry of the 
messengers, together with gauge invariance, would allow only for the extra coupling 
Q 0 i( 5 ,d) 02 ,( 5 ,l )5 which is forbidden by the discrete R-symmetry in Table 1. 
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The boundary conditions induced at the messenger scale by this SUSY breaking 
sector includes gauge mediated contributions to the scalar masses and the contribu¬ 
tion arising because of the new interaction A. The soft masses for this model are 
obtained in the Appendix A, here we report simply the results, dehning A = ^. The 
gauge mediation contribution is 



2 

(167r^)2 


YU9‘(2AV.(T)) 


(2.4) 


where fs{x) is the usual minimal gauge mediation function for sfermions (see e.g. [52]) 
fs{x) = 1 + ^ + 0{x^). The contribution to the soft masses induced by the A 
deformation are 


m 


2 _ 
UiUj - 


1 

2567r4 


Cij A djj 



= 


mrr = 
tiu 


dij}? 

2567r4 

Sdi/X^ 

2567r4 


iy^.c.y)ij A^ 
Triy^.c.y) A^ 


2 Crg^r 


A^ 


du ,2t / ^ 


(2.5) 


where we dehned the matrix Cij = CiCj, du = 2 and d^ = 4, Ci = 2/5, C 2 = 0, Us = 4. 
We give the complete expression for the one loop function h{x) in Appendix A. It 
can be approximated by h{x) = 1 -|- -f- 0{x'^). 

The up type right squark gets off diagonal contributions whose flavour structure 
is determined by the matrix c^, consisting of a two loop contribution (which can 
be positive or negative depending on the value of A) and a one-loop negative con¬ 
tribution, which can be relevant for not too small values of Besides the up-type 
squark, the rest of the soft masses contributions are determined by the Yukawa cou¬ 
plings yij. In particular the contribution to rug.^., that would have implied otherwise 
strong bounds from flavor observables, is projected along the Yukawa couplings. 

Moreover, the fact that the hidden sector respects a discrete R symmetry implies 
that we have not generated any A-term. This is a welcome feature in perspective 
of possible large off-diagonal contributions induced by the A deformation. Indeed, 
large off diagonal A-terms would have raised two problematic issues. First, large off 
diagonal A-term involving the hrst generation can lead to large neutron EDM, above 
the current experimental bound, de facto excluding the model, as we previously men¬ 
tioned. Second, large off diagonal A-terms, together with gaugino masses, induces 
corrections to fermion masses which can be incompatible with the actual value of the 
light quark masses. On the other hand, the absence of A-terms implies that in order 
to get the correct Higgs mass we will have to consider at least one stop to be quite 
heavy. 

In order to induce non vanishing gaugino masses and sizable scalar masses we 
consider also in addition to the previous SUSY breaking sector another supersym- 
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metry breaking sector inducing general gauge mediation contribution, with different 
gaugino and scalar SUSY breaking scales, respectively and As 





9r 


Ibvr^ 
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A 


(IGtt 


2\2 


G 




2 

s- 


( 2 . 6 ) 


For simplicity we assume that such contribution is induced at the same messenger 
scale M as above, which hence sets the range of scales of the renormalisation group 
(RG) flow. 

In summary, the total contribution to the soft masses in the complete model is 
given by adding expressions (2.4, 2. 5, 2. 6). 


We explored the parameter space of this model and the implication for flavor 
observables. We implemented the model in SARAH [53, 54] and generated the spec¬ 
trum and compute the flavor observables using SPheno [55]. We also computed the 
contribution to Dq mixing, which is sensitive to the mixing between the first and 
second generation squarks, and the neutron EDM, using the SUSYFLAVOR code 
|56]. 

As a large contribution to Dq mixing arises through hadronic long-distance 
physics [58, 59], which is plagued by large theoretical uncertainties, for the limit 
setting we only impose an upper bound, i.e. its measured central value, on the short- 
distance SUSY contributions. The flavour observables that we checked, together with 
the bounds that we applied, were taken from [57] and are summarized in Table 2. 

In order to establish general low energy physics constraints on the model, we 
performed a scan by fixing the values of A^, As and A and varying the messenger 
mass M, the deformation size A, and the flavour direction c of the deformation. In 
the left plot of Fig. 1 we show the result in the (A, A/M) plane. The points respecting 
flavour constraints are shown as circular dots, while the crosses are points violating 
flavour observables. The red points are scenarios where the lightest up type squark is 
lighter than 1.5 TeV. The Higgs mass is correct (within the errors) on all the points 
shown in the plot, due to the large value of A^. Since we are marginalizing over 
the flavour direction, there are overlapping points which have the same size of the 
deformation but which can or cannot satisfy the flavour bounds. However the plot 
is useful in order to understand the effect of the deformation on the soft spectrum, 
independently on the direction c, as we explain now. 

The main effect of the deformation is on the mass squared (see eqn. (2.5)), 
which can be negative (or very small) at the messenger scale and hence positive 
but small at the EW scale, in the region of moderate A and large ratio A/M. In 
particular the final spectrum is very sensitive to the value of the ratio A/M, because 
it determines how large is the one loop negative contribution to the soft mass mf;.. 




Flavour Observable 


Imposed limit 


Source 


Tool 


BR{B—>-Xs^)sm 

B R(B'^—>- fi/i) g M 

BR(B°^tifi)sM 

(AMs, )sM 

_ {^Mb^)sM 

BR{D^fii') 

BR{D—^iii/)sm 

BR(Ds^^iu) 

BR(Da^^iv)sM 

BR{K^Im^ 

BR{K—>-iii/)sm 

BR{B^s^i^i) 

BR{B—>-snii)sM 

BR{B^K fifj.) 

BR{B^Kij,ij,)sm 

_ (^k)sm _ 

AMk 

_ (^Mk)sM 

AMd 

Neutron EDM \d 


0.84, 1.16] 

[0.87, 1.08] 

[0.90, 1.06] 

[0.90, 1.17] 

[0.85, 1.13] 

[0.95, 1.01] 

[0.95, 1.03] 

[0.99, 1.01] 

[0.12, 1.87] 
[0.875, 1.125] 
[0.68, 1.34] 
[0.997, 1.003] 

< 8.82 X 10-^5 GeV 

< 2.9 X 10“^® (e cm) 


HFAG [60] 


HFAG [60] 


HFAG [60] 


HFAG [60] 


HFAG [60] 


HFAG [60] 


HFAG [60] 


HFAG [60] 
HFAG [60] 
HFAG [60] 


[61 


|61| 


|62| 


|63| 


SPheno 

SPheno 

SPheno 

SPheno 

SPheno 

SPheno 

SPheno 

SPheno 

SPheno 

SPheno 

SPheno 

SPheno 

SUSYFLAVOR 

SUSYFLAVOR 


Table 2: Flavour bounds imposed during the scan on the model parameter space. 

at the messenger scale. The allowed region in Fig. 1 is determined by the fact that 
for large values of A/M the negative one loop contribution to the up-type squark 
mass (independently on the flavour direction) renders tachionic the squark eigenstate 
aligned with c, and the spectrum is rejected. Indeed, the points where there is a single 
light right-handed up squark with mass smaller than 1.5 TeV, the red points, are at 
the border of the allowed region. The two loop contribution to the up squark mass in 
(2.5) can be positive or negative, depending on the value of A. This determines the 
shape of the allowed region, and the location of the throat, which is approximately 
where the two loop correction in (2.5) changes sign. 

Note that in the red region only one mass eigenstate ui is light and much lighter 
than all the other squarks. This mass eigenstate will be aligned along the vector c 
in flavor space, as the deformation A (if we neglect the effects of the GKM mixing). 
This can be easily understood by observing that the gauge mediation contributions 
are diagonal in flavor space while the A deformation contributions are proportional to 
the matrix Cij = CiCj, which has by construction c as eigenvector. Hence the lightest 
right up-type squark will be a mixture of right handed stop, scharm and up-squark 
in the combination ^ 


ui — Ui^ur + U 15 CR + Uietn with |Hi4|^-|-— 1 (2.7) 

^We used here the SLHA2 notations [64] (i.e. 4 <->■ ur ,5 ^ cr ,6 ^ tR). 
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Figure 1: Left: results of the scan on the parameter space of the model on the (A, A/M) plane. 
We fixed Ag = 1 x 10® GeV, Ac = 1.9 x 10® GeV (corresponding to gluino mass around 1.5 TeV), 
A = 5 X 10® GeV. We scanned over A, ci,C 2 ,C 3 and M. The red points have one squark with 
mass smaller than 1.5 TeV. Right: we focus on a region with light squark: we fixed A = 1 and we 
restricted M to be such that A/M ~ 0.2, varying freely on (ci, C 2 , C 3 ); we show only points having 
one light squark with mass smaller than 1.2 TeV, and we plot in the physical mixing angle plane 
(|17 i4|, IG 15 I), where 17i4 is the up component and Ui^ is the charm component of the lightest squark 
(see eqn (2.7)). All points displayed have mh = 125 ± 2.5 GeV. 


where \Ui^\ ~ ci, |t/i5| ~ C2, \Uiq\ ~ C3. 

The contribution to the other squark mass rug is also negative, but it is a two loop 
effect and it is compensated by the large value of the gauge mediation contribution 
since we have considered large A 5 . This is also a welcome feature to alleviate possible 
sources of flavour violation in the Q sector descending from the deformation A, since 
the gauge mediation contributions are flavour diagonal. The other generic effect of 
the deformation is that is very small at the scale M, and hence it becomes very 
large and negative at the EW scale. As a consequence, since the EWSB condition 
sets approximately p ~ —the /i parameter is very large at the EW scale. Hence 
the Higgsino is typically quite heavy. 

In the plot on the right of Fig. 1 we focus on points with one light up-type 
squark with a mass smaller than 1.2 TeV. We fix A = 1 and set A/M in the range 
around 0.2. Eventually we vary over c. The plot shows the result of the scan on the 
\ Ui 5 \ mixing angle plane (remember that Yhi Uu = !)• 

The red crosses are points which are not allowed because they do not respect 
flavour observables, while the circular dots are viable points. As mentioned above, 
the main flavour violating effect induced by the deformation A is in the RR up 
squark mass matrix, precisely the one loop negative contribution in (2.5). There is 
no off diagonal contribution to the A-terms, hence no LR mixing, and the two loop 
contributions to the LL squark masses induced by the deformation A are negligible 
compared to the flavour diagonal gauge mediated one ( 2 . 6 ) for all the scanned points. 
As a consequence the only relevant flavour constraint in Fig. 1 is actually coming 
from Dq — Dq mixing. The plot shows that off diagonal contributions to RR up 
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Figure 2: Higgs mass constrains on the supersymmetric spectrum, assuming parameter as in Fig. 
1, (i.e. As = 1 X 10® GeV, Aq = 1.9 x 10® GeV, A = 5 x 10® GeV), and where we furthermore set 
A = 1 and C 2 = 0. The red region has mt = 125 ± 2.5 GeV, the yellow region mu = 120 ± 2.5 GeV, 
the black region even smaller mh- The mass of the lightest squark mass is denoted with rriu-^ and 
? 7 i 4 is the sup-stop mixing angle (here [/15 ~ 0 ). 


squark masses can be compatible with low energy flavour contraints if the mixing 
between the first and the second generation is not large. Hence in this model we 
can realize scenarios where the lightest squark is a highly mixed state in the flavor 
basis, either sup-stop or scharm-stop (essentially the circular red points along the 
two axes). 

In the plots, all points have a viable Higgs mass, which is a consequence of the 
fact that at least one of the two stops is very heavy, with a mass set mainly by A 5 . 
However the Higgs mass constraint implies that the deformation cannot lower too 
much the mass of the lightest eigenstate, if this has a significant component in the 
direction 3 in flavor space. Hence very low mass values for the lightest squark are 
allowed only if they are not aligned along the third family. This explains the hole 
of points in the region of small Ci and C 2 , where the deformation is mainly aligned 
along the direction 3 in flavour (note, we have fixed A = 1 and that c\ + 0 ^ + c\ = 1). 
For a recent treatment of the Higgs mass formula in a fully flavored MSSM see e.g. 
165], 

This effect can be qualitatively estimated by considering Fig. 2, where we fix 
C 2 = 0 and we plot the contours for the value of the Higgs mass (now allowed to take 
also small values) as a function of the lightest squark mass and of the mixing angle. 
The region with viable Higgs mass is denoted in red. The lighter the squark is, the 
more it should be aligned along the direction 1 in flavour space, and hence less in 
the direction 3 in flavour, in order to have a large enough Higgs mass. 
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1.9 X GeV 


M 

2.31 X 10® GeV 


A Cl 
1.0 0.75 


rriu^ 

\Uia\ 

\Uw\ 

\Ui,\ 


446 GeV 

0.706 

~ 0 

0.707 

m-g 


rrig 

rrih 

m3i2 

1.56 TeV 

526 GeV 

260 GeV 

124 GeV 

1/k X 612 eV 


A 


G 


2.5 X 10*^ GeV 


M 


2.315 X 10^ GeV 


A 


1.0 


Cl 


0.75 


rriu^ 

\Ul4\ 

\Ul5\ 

\Uw\ 


657 GeV 

0.706 

~ 0 

0.707 

m-g 



rrih 

m3i2 

1.99 TeV 

688 GeV 

345 GeV 

124 GeV 

1/k X 614 eV 


Table 3: Two example of viable spectra. The supersymmetry breaking scales for the scalar sector 
are fixed to A 5 = 1 x 10® GeV and A = 5 x 10® GeV. The parameter C 2 is set to zero, hence 
c| = 1 — cf. Goncerning the rest of the spectrum, for both benchmarks: the sleptons are at least 
heavier than 1 TeV; the other squarks and the higgsinos are further heavier, with masses larger 
than few TeV’s. 

From Fig. 2 it is evident that we can easily get viable points with a correct 
Higgs mass, satisfying flavor constraints, and having a very light and highly mixed 
right-handed up squark. In the scatter plots we hxed to a representative value, at 
the border of the LHG8 exclusion reach for the gluino mass; however, the qualitative 
feature of the typical spectrum does not change by raising further the gaugino mass 
scale. In Table 3 we show two prototypical benchmark points, with different gluino 
mass, both presenting one lightest squark eigenstate, maximally mixed between stop 
and sup gauge eigenstates. 

Summarizing, the deformation presented in this section provides a model real¬ 
ization of scenarios with a light squark with a large stop-sup or stop-scharm mixing. 
The drawback of this model is that the Higgs mass is obtained by requiring the 
other squarks in the spectrum to be heavy, implying a certain unavoidable amount 
of tuning in order to have one light squark state. 

2.3 Mass Spectrum 

In this section we discuss the typical spectra that are generated by the model pre¬ 
sented in the previous section. 
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We are interested in regions of the parameter space where there is a very light 
right-handed up squark, which is highly mixed in flavor. In the following we focus on 
the case of large sup/stop mixing, assuming a small scharm component, consistently 
with the points shown in Fig. 2. We have shown in the previous section that 
the complementary case of scharm/stop mixing can be also easily generated. The 
phenomenology of such case has been recently discussed in [ 6 ]. 

The supersymmetry breaking scale determining the gauge mediation contribu¬ 
tion associated with the scalars (A 5 ) is typically large in order to have the other 
squarks in the few TeV range. As a consequence, also the sleptons are very heavy 
and essentially decoupled from collider physics. The higgsinos are also heavy and 
decoupled, following the argument explained in section 2 . 2 . 

Hence the gauginos and the gravitino are the only other supersymmetric particles 
that can play a role in the collider phenomenology of these models. In the model 
formulation, we made the assumption that we have only one gaugino mass scale A^. 
This implies that in the benchmark points shown in Table 3 the gaugino masses 
respect the following relation involving the gauge couplings M 3 : M 2 : Mi = : 

g 2 : gf. However, we expect that this relation can be relaxed considering a complete 
GGM parameter space [ 66 , 67], effectively disentangling the Bino from the Wino and 
from the gluino mass. Note that modifying the gaugino mass at the messenger scale 
will affect the running of the sfermion masses, and hence could modify quantitatively 
but not qualitatively the results we obtained above. 

The LSP is always the gravitino, whose mass is given by 


m3/2 


Mv/A2 + Max(A^,A|) 


i/SkMpi 

and all the sparticles have universal decay to it with the formula 


r(x ^ XG) ~ 


m 


X 




( 2 , 8 ) 


Here A; is a factor which is smaller than 1, taking into account that the fact that 
the supersymmetry breaking parameters coupling to the messengers can be smaller 
than the supersymmetry breaking scale of the complete model. For instance, if the 
supersymmetry breaking scales A, Ac, As are generated at one loop in some model 
of dynamical supersymmetry breaking, the factor fc is a one loop suppression. In the 
Table 3 we show also the value of the gravitino mass for the prototypical benchmark 
point. Depending on the value of fc, the decay of the NLSP can be displaced or 
longlived. In the following we assume the second case, such that the Bino decay does 
not play any role in collider signals. 
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ui + r/w* 


Ui = + in) 


t/u 


Xo 


Figure 3: Maximally Mixed Sup-Stop (MMUT) scenario: the simplified model considered in the 
collider analysis. 


3 Simplified model 


Based on the previous analysis of the parameter space of the model, we here define 
the simplihed spectra that we investigate in the phenomenological study of the rest 
of the paper, and we set our benchmark points. The only light supersymmetric states 
are a highly mixed right handed up-type squark, that we denote with -ui, with masses 
varying between 400 and 1000 GeV, and a pure Bino neutralino Xoi with mass varying 
between 200 and — 200). The gluino, which is relevant since it participates in 
the production mechanism of the light squark, and can also be directly produced, is 
hxed to 2 TeV. We focus on the case with maximal stop-sup mixing, i.e. the light 
squark state is a perfect mixture of sup and stop right-handed squarks 


— -^{ur -|- Ir). 


The simplihed model, that we dub as Maximally Mixed Sup-Stop (MMUT) scenario, 
is shown in Fig. 3, and is motivated by the spectrum structure that we described in 
the previous sections. In general, the Wino is considered to be heavy enough not to 
participate either in the squark production processes or in possible squark decays. ^ 


3.1 Branching Ratios 

We are interested in production of the lightest squark eigenstate and the resulting 
decay chain. Hence here we study the possible decay channel for ui. 

The decay of Ui to the gravitino (given also by the universal formula (2.8)) is 
always suppressed compared to decays via gauge couplings. If the mass difference 
between ui and xo is smaller than the top mass, ui will decay to m -|- xo with a 
branching ratio which is almost 1, as soon as there is a m component inside ui. In 

^Actually EW production via Wino and also Bino exchange can increase the cross section at 
LHC of a few % level, but we do not consider it in the following. 
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this case its signature will resemble the one of a light up-type squark. This can lead 
to interesting bounds in some regions of the squark-neutralino plane [7]. 

A more unusual scenario can be realized when the mass difference between the 
Ml and the xo is larger than the top mass‘d. We focus on this parameter region in the 
following sections. Here the two competing decays are ui —)• uxo and ui —)■ txo- 

The branching ratio for decays into top quarks is slightly suppressed even in 
the case of maximal mixing due to the phase space suppression from the large top 
mass. The formula for the branching ratio can be extracted by the analytic results 
of [68, 69], and in the case of maximal mixing reads 


BR[ui txo] = 


m2o, m2](m|i - - mj) + (m|^ - J 


2 ’ 


(3.1) 


where 


-2m? 


mt 


xo 


2 m? ^m: 


2m2^mj. 


In Fig. 4 we plot the branching ratio of Eq.(3.1) on the ('Ui,xo) mass plane. Note that 
the branching ratio into top is still sizable even in region of moderate compression of 
the spectrum. For instance for rriu^ = 450 GeV and = 200 GeV the branching 
ratio into top is still between 30% and 40%. For an even more compressed region 
the branching ratio into tops drops abruptly. 

Given the spectrum of the simplified MMUT scenario in Figure 3, at the LHG 
also the associated production of gluino and the lightest squark state can play a 
significant role. The gluino can decay either as g ^ UiU* or as ^ Uit* (and the 
conjugate processes), hence giving extra jets or extra tops in the final state. Since 
the mass difference between the gluino (at 2 TeV) and the lightest squark ui is large 
(at least larger than 1 TeV on our parameter space), the phase space suppression 
induced by the top mass is always negligible. Indeed we checked that the branching 
ratio of the gluino decay into squark and up or top quarks is always proportional to 
the mixing angle. Hence in the MMUT model of Figure 3 the gluino will decay 50% 
to up-quark and wi, and 50% to top-quark and ui. 


3.2 Benchmark points 

We define five benchmark points on which we concentrate in the collider analysis, all 
with gluino mass fixed at 2 TeV, and with different squark and neutralino masses. 
They are reported in Table 4. The first one has very light squark but quite com¬ 
pressed spectra, with m^j — m^g = 250 GeV. The others have larger squark masses 
and light or moderately heavy neutralino. In the next sections we will discuss the 
exclusion limits on such benchmark at LHG8 and their distinctive signatures for 
LHG14, characterized by tops in the hnal state. 

^Another interesting possibility is if xo is heavier than ui, which can then decay only to the 
gravitino. 
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BR[mi^ t^o] 



Mn, 

Figure 4: Branching ratio BR(ui —>■ txo) as a function of the ui mass and of the neutralino mass. 



m-g 


^XO 

Ui mixing angles 

Benchmark Point 1 

2 TeV 

450 GeV 

200 GeV 

|f/l4 

= |f/i5 

1 

“ ^/2 

Benchmark Point 2 

2 TeV 

700 GeV 

200 GeV 

|Gi4 

= \Ul5 

_ 1 
\/2 

Benchmark Point 3 

2 TeV 

700 GeV 

400 GeV 

|f/l4 

= |f/l5 

_ 1 

U2 

Benchmark Point 4 

2 TeV 

950 GeV 

200 GeV 

|f/l4 

= |f/i5 

1 

“ P2 

Benchmark Point 5 

2 TeV 

950 GeV 

400 GeV 

|Ul4 

= |Gi5 

1 


Table 4: Benchmark points considered in the collider study of the simplified MMUT model. 


4 Collider signatures 

4.1 Production modes and cross sections 

Given the simplified spectrum of the MMUT scenario in Fig. 3, at the LHC we 
expect the following production modes: 

pp —>■ uiu\ , pp —)• 'Ui'Ui , pp —)• u\g (4.1) 

Cross sections for squark-antisquark, squark-squark and gluino-squark production at 
LHC8 and LHC14 are shown in Fig. 5, computed at LO using MadGraph 5 [70]. 

Note that gluino pair production, gluino-antisquark production and antisquark 
pair production do not contribute significantly to the SUSY production modes, be¬ 
cause of the PDF suppressions. To give an estimate of these contribution and com¬ 
pare with the cross sections depicted in Figure 5 we evaluate these cross section with 
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Figure 5: Cross sections in the MMUT model for squark-antisquark, squark-squark and gluino- 
squark production at LHC8 and LHC14 as a function of the squark mass. The gluino mass is fixed 


to 2 TeV. 


rriu = 400 GeV and rrig = 2 TeV. The cross section for gluino pair production is 
a ~ 10“^ fb and a 0.7 fb at 8 TeV and 14 TeV respectively, while the cross section 
for gluino anti-squark production is a ~ 0.02 fb and cr ~ 1 fb. For the anti-squark 
pair production cross section we hnd a ~ 0.2 fb and cr ~ 1.6 fb at 8 TeV and 14 TeV 
respectively. Hence they are negligible and we will not consider them in our analysis. 

We stress that the production modes of squark-squark and gluino-squark are 
sizeable only because of the up component in the lightest squark, and thus they 
are weighted by the mixing angle. For instance the squark-squark production is 
proportional to | Vi 4 |^, which is 1/4 in our maximally mixed scenario. This production 
mode vanishes in the scenario where the lightest squark is purely stop-like, and is 
maximized when the lightest squark is a pure up squark. 

Including the possible decay modes, the relevant production mechanisms are 
depicted in Figure 6. Depending on the decay of the gluino and of the lightest squark, 
the hnal state can be composed of up quarks or of top quarks and neutralinos. Since 
we consider the neutralino to be stable in terms of collider time scales, it will give 
rise to missing energy signatures. 

One can generate the following signals from squark-antisquark pair production 
(right Fig. 6) 

pp uiul -)■ {jjxoXo, jtXoXo, tt*XoXo) (4-2) 

the following hnal states from squark-squark pair production (left Fig. 6) 

pp -)■ uiui -)■ (jjxoXo,, jtXoXo, ttxoXo) (4.3) 

and the following signatures from squark-gluino associated production (bottom Fig. 

6 ) 


pp uig {jjjXoXo, jjtXoXo, jjt*XoXo, jtt*XoXo, jttxoXo, ttt*xoXo) (4-4) 
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Figure 6: Top left: squark pair production (note that this process is produced through t- 

channel exchange of gluino). Top right: squark antisquark pair production. Bottom: squark gluino 
production (also the conjugate decay process for the gluino decay should be included). 


Besides the usual supersymmetric signatures of jets plus Ifij, or of top-antitop pair 
plus l^rp, among the possible final state we can find the single top (hence single 
lepton) and same sign tops (hence two positive same sign leptons). Especially this 
second signature, i.e. same sign tops, is a unique consequence of the maximal flavour 
mixing in our MMUT scenario, and can be considered as a probe of the mixing angle 
of the lightest squark state. 

4.2 Constraints from LHC8 

The supersymmetric processes described above contribute to final states with jets 
and missing transverse energy, and are possibly probed by the jets plus Ipj, searches 
of ATLAS [71] and CMS [72] at 8 TeV. Analogously, presenting top pair and Ipj, in 
the final states, they could be also probed by the standard stop searches of ATLAS 
[73] and CMS [74], 

We thus have to verify that our benchmark points in Table 4 are not excluded by 
existing LHC8 searches. In order to ascertain the viability of our benchmark points 
we produce samples with MadGraph 5, shower them with Herwig++ [75], and process 
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them using CheckMate [76], using all the available LHC searches, including jets plus 
Iprp [71], searches exploting the ut variables [77], and single-lepton stop searches [73]. 
We hnd that all hve benchmark points in Table 4 are still not excluded by LHC8 
searches. The benchmark points which are at the border of the exclusion are the hrst 
and the second one in Table 4. However, with an overall K factor as large as® 1.7, 
the CheckMate analysis concludes that these two benchmark points are still allowed. 
The stronger constraints come from the ATLAS jets plus Ipj, search [71]. 

This result could be surprising given that the SUSY cross section is enlarged 
by the contribution of the squark-squark production mode with respect to the un¬ 
mixed squark scenario, where the only relevant contribution is coming from squark- 
antisquark production (the gluino-squark associated production is negligible at 8 
TeV in the small squark mass region). However, there are two effects which both 
contribute in making the signal of the MMUT model difficult to exclude. First, the 
cross section enhancement due to the squark-squark channel is only mild, given that 
this production mode is suppressed by the mixing angle to the fourth power; this 
results in a factor of 1/4 suppression, as mentioned above (see Figure 5). Second, 
the different decay modes of the light squark eigenstate, including top quarks, make 
the signal less clean than pure jets plus topologies. This is particular effective 
in the hrst benchmark point (with mu^ = 450 GeV and = 200 GeV). In this 
case the branching ratio of the lightest squark into top is still sizeable (see Figure 
4), but the spectrum is moderately compressed and the resulting tops will be quite 
soft, effectively reducing the efficiencies. 

In order to give an intuitive explanation of these features, in Figure 7 we plot 
the total cross section of our supersymmetric production modes, weighted with K- 
factors taken from [78, 79], and the cross section exclusion limit extracted from GMS 
[74] in the case of neutralino mass at 200 GeV. This plot is an oversimplihcation, 
since it does not take into account the different decay modes and efficiencies in the 
various channels. Nevertheless it allows to get some insight on the LHG8 reach for the 
MMUT simplihed model and how the mixing angle plays a crucial role in reducing the 
experimental bounds. We observe that the region of squark masses between 500 GeV 
and 680 GeV is excluded, with the 450 GeV and the 700 GeV benchmarks just above 
the experimental reach, consistent with our results from GheckMate. To provide a 
further quantitative estimate of the effectiveness of the mixing angle suppression, 
in the plot we also show the total production cross section for a right handed up 
squark, not mixed. In this case the same sign production mode is not suppressed by 
the mixing angle and the total cross section is considerably larger, leading potentially 
to very strong bounds. In our plot, the exclusion limit for scenarios with a pure up- 
type right handed squark seems to reach very large values for the squark mass; this 

® Considering the relevant production modes involved in our process, a K-factor of 1.7 is the 
largest possible value at 8 TeV [78, 79]. 
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Figure 7: Simplified exclusion plots based only on the SUSY production cross section. The 
neutralino mass is fixed to — 200 GeV and the gluino mass to 2 TeV. The dashed, blue 
curve is the LO total SUSY production cross section for the MMUT model, including all relevant 
production channels (miu*, wifti, uig). The solid, blue curve is the same total cross section for the 
MMUT model, but re-weighted by the appropriate K-factors; we used = 1-6, = 1-1 

and Ku^g = 1.7, taken from [78, 79]. The solid, dark red curve is the exclusion cross section 
extracted from the CMS analysis [74] for neutralino mass at 200 GeV. The dotted, red curve, for 
comparison purposes, is the cross section for the case of a pure (right handed) up squark, also 
re-weigthed with the appropriate iiT-factors. 


is a consequence of the oversimplification of this estimate, for a complete recasting 
and evaluation of the LHC reach see [80]. 

Other searches that could potentially constraint the model are the standard 
supersymmetry searches for stops [73, 74]. However, note that the total production 
cross section in our scenario is smaller than in the case of a pure light stop. Indeed, 
even if in addition to the squark-antisquark production mode we have also the squark- 
squark production mode, the entire cross section should be multiplied by a factor 
of BR{ui —>■ txo), in order to require at least one top in the final state. One can 
easily check that for instance on the first benchmark point (with = 450 GeV 
and = 200 GeV) this implies that the cross section times the branching ratio 
for a single top in the MMUT model would be eventually lower than the single 
top production in a model with a light pure (right handed) stop with the same 
mass. Hence, these simplihed estimates already provide an understanding which is 
consistent with the robust result obtained through Gheckmate. For squark masses 
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larger than 700 GeV we do not expect existing stop searches to be relevant. 

For sqnark masses smaller than 450 GeV, the branching ratio into tops drops 
signihcantly, to get completely negligible already at a sqnark mass of 420 GeV (see 
Fignre 4). Moreover, the spectrnm is considerably compressed and hence weakly 
constrained by LHG searches. Hence in the MMUT model there is a small window 
for sqnark mass aronnd 450 GeV which is still allowed by LHG8 searches and which 
present tops in the hnal state, thanks to the pecnliar properties associated with 
maximal flavonr mixing. In the following, we hnd therefore interesting to inclnde 
this blind-spot benchmark point in onr analysis and discnss its phenomenology at 
LHG14. 

The above discnssion is valid for a nentralino mass of 200 GeV. For lower nen- 
tralino masses, the bonnd on the sqnark mass is stronger. Gonversely, a heavier 
nentralino mass leads to weaker bonnds. 

To snmmarize, the hve benchmark points for the MMUT model we dehned in 
Table 4 are still allowed by actnal LHG8 searches. Nevertheless, we expect that in 
the next LHG rnn at higher energy the standard SUSY searches for light sqnarks 
or stops conld eventnally probe them. In the following we concentrate on the new 
signatnres associated with large flavonr mixing, and access the LHG14 reach for snch 
signatnres. These are new interesting channels to look for snpersymmetric scenarios, 
and at the same time represent a powerfnl strategy to test the mixing property of 
the light sqnark. 

The collider signatnres distinctive of onr simplihed models are: 

1. Single-top signal arising from all prodnction modes (4. 2,4.3,4. 4) 

2. The new signatnre of the snp-stop mixing, i.e. same sign-tops, arising from the 
prodnction modes (4.3) and (4.4) 

In the next section we will stndy the prospects for these two new channels at LHG14 
for the hve benchmark points of the MMUT simplihed model. 

5 Top signatures at LHC14 

5.1 Event Generation and Reconstruction 

We generate all events nsing leading order MadGraph 5 [70] with the NNPDF 2.3 [85] 
set for the parton distribntion fnnctions. Upon hard level process generation, we 
fnrther shower the events nsing Pythia 6 [86]. We match the backgronnd event 
samples, where relevant, to extra jets nsing the MLM matching scheme [87], with 
QCUT=30 GeV and xqcut = 20 GeV in case of top prodnction and QCUT=15 GeV and 
xqcut = 10 GeV for prodnction of weak bosons. Onr analysis inclndes detector 
ehects on event reconstrnction, where we ntilise Delphes 3 [88] with the defanlt 
GMS settings for event reconstrnction, Ytagging efficiencies and lepton isolation. 
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In the following sections we consider signatures of our model in the single top (i.e. 
I + + b) channel as well as the same sign positive top (i.e. channel, 

requiring us to consider a range of SM background channels. For the purpose of 
studying the collider reach in the single top channel, we consider tt and W + bb 
events matched to one extra jet as well as tW events, where we require at least one 
lepton of unspecihed charge at hard process generation level. 

Our SM background samples for the same sign lepton analysis consist of ti and 
W + bb events matched to one extra jet, where we require one positively charged 
lepton at generator level. In addition, we also consider rare SM processes where we 
include production of tiW,tiZ, ZZ,W^W~^W~ and WZ. Here we require at least 
one positive lepton in the hnal state. 

In order to improve our estimates of background event yields, we normalise the 
ti production cross sections to the NNLO+NNLL value of Ref. [89], while we assume 
a conservative K-factor of 1.4 for W + bh and single top production and 1.3 for rare 
SM processes. For the purpose of signal generation, we always assume a K-factor of 

1.1 for uu, 1.4 for ug and 1.5 for uu* production [78, 79]. 

5.2 Single-top Channel 

Our analysis of the single top channel {i.e. I + Ifij. -|- b) is inspired by previous work of 
Ref. [82]. We include modihcations to optimise the analysis for the high luminosity 
LHC. The single top event topology of signal events results from all the processes 
depicted in Figure 6. The pair production of either same or opposite sign u contribute 
to this hnal state when one of the squarks decays into a xo and a w-jet, and the other 
decays into xo and a t quark®. Analogously, the gluino squark associated production 
generates single top topology when either the gluino or the squark present at least 
one top in its decay chain. 

Fig. 8 shows examples of interesting kinematic distributions for signal and back¬ 
ground events. Signal events in the single top channel are characterized by large 
missing energy compared to SM backgrounds, as well as large transverse mass^. The 
transverse mass distributions of the SM backgrounds display a suppression around 
the W mass, as the only source of signihcant missing energy and hard leptons is the 
decay of the W boson. Conversely, a signihcant contribution to missing energy in 
the signal events comes from xo, allowing the transverse mass distribution to extend 
to much larger values. 

The degree of squark-neutralino mass degeneracy has a large ehect on the shape 
of transverse mass and missing energy distributions. Lower panels of Fig. 8 il- 

®We note that mixing between sup and stop can give rise to loop-induced direct top 
production[10]. This contribution could potential increase our signal rate, however, after apply¬ 
ing the analysis cuts of Eq. (5.1) direct top production is irrelevant. 

^Here we define rriT = 
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Figure 8: Kinematic distributions of signal and background events in the I + fj, + b channel. All 
distributions are normalized to unit area and assume no cuts. The signal distribution in the upper 
panels shows a benchmark point of rriu = 450 GeV and = 200 GeV. In the bottom panels, the 
values in the labels represent (m^, 


lustrate the point for the benchmark mass of = 700 GeV. The benchmark 
point with = 400 GeV, leads to much softer Ipj, and rriT spectra compared to 
= 200 GeV in uu and uu* production, while the effect is much milder in the ug 
production. The effect suggests that a cut on or Ipj, which would be appropriate 
to isolate the signal in the compressed mass region, might not be optimal in the 
non-compressed spectrum. For the purpose of illustration, here we will focus only 
on cut selection criteria which can better probe the uncompressed mass spectrum 
scenario. 

In order to improve the purity of the event sample, we impose the following set 
of kinematic cuts: 

Ni{pt > 20 GeV, g < 2.5) = 1, Ni,{pt > 20 GeV, p < 2.5) = 1, 

> 250 GeV, rriT > , (5.1) 

where h and I refer to 6-tagged jets and isolated leptons respectively, and we choose 
= (210, 310,410) GeV for mu = 450, 700, 950 GeV respectively. 
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rriu = 450 GeV, = 200 GeV, rrig = 2 TeV 


1 l/jrp -\- b 

uu 

ug 

uu* 


tt 

W + bb 

W + t 

S/B 

S/y/B (300 fb-^) 

Ni = l,Nh = l 

2.6 

39.0 

1.1 

4.1 X 10^ 

420.0 

405.0 

1.0 X 10-3 


3.6 

> 250 GeV 

0.67 

0.82 

8.3 


238.0 

16.0 

9.3 


0.037 


10.0 

mr >210 GeV 

0.27 

0.40 

3.4 


12.0 

0.15 

< 0.1 


0.32 


20.0 



rriu = 

950 GeV, 

= 400 GeV, mg 


2 TeV 



1 1/jrp -\- b 

uu 

ug 

uu 

* 

tt 

W + bb W + t 

S/B 

S/^/B (300 fb"^) 

Ni = l,N, = l 

0.37 

0.36 

0.48 

4.1 X 10 

^ 420.0 

405.0 

2.9 X 10“ 

5 

0.10 

> 250 GeV 

0.26 

0.30 

0.34 

238.0 

16.0 

9.3 


3.4 X 10“ 

3 

0.96 

tut > 410 GeV 

0.093 

0.12 

0.13 

3.1 

< 0.1 

< 0.1 

0.11 


3.4 


Table 5: Example cutflow in the I + f!rp + b channel. The entries show cross sections in fb after 
each consecutive cut. hE+jets and tW channels includes a generation level cut of Ipj. > 80 GeV in 
order to improve the statistics in the signal region. 


( GeV) 

' 450 

700 

950 

200 

400 

' (0.32, 20.0) 

(0.24, 11.0) 
(0.11, 5.3) 

(0.12,3.8) 

(0.11,3.4) 


Table 6: Summary of reach for benchmark points in the single top channel. The table entries 
show S/B and S/'/B at 300 fb“^ respectively which can be achieved at LHC Run II in the single 
top channel. The masses of squarks and neutralinos are listed on in the topmost row and leftmost 
column respectively. All results assume rrig = 2 TeV. 


Table 5 shows an example cutflow for a benchmark signal points of = 450,950 
GeV and = 200,400 GeV. Requiring exactly one lepton and at least one 6-tagged 
jet is sufficient to bring the signal to background ratio {S/B) to levels of ~ 10“^ for 
lighter squark masses, but only to ~ 10“^ if the squark mass is ~ 1 TeV. A cut on 
missing energy results in a factor of ~ 100 improvement in S/B, while the additional 
cut on rriT improves S/B by an additional factor of 10 at a 50% signal loss. Our 
results, summarised in Table 6, show that the mu = 450 GeV, m^^ = 200 GeV 
benchmark point is discoverable with the signal significance of S/\fB ^ 5 with 
L = 300 fb“^ of integrated luminosity at LHG Run II, while with the same amount of 
data we should be able to rule out our model for squark masses of roughly < 1 TeV 
A potential discovery of a signal in the I + J/Jj, + b channel would give indirect 
evidence for the existence of supersymmetry, but would not provide information on 
the degree (if any) of the u — t mixing, as even then minimal flavor-conserving SUSY 
models predict signals in the single-lepton channel. Measuring additional channels 

®The cuts we chose for this analysis are somewhat optimised for high luminosity LHC. Relaxing 
the rriT and cuts at lower integrated luminosities is also likely to be more efficient in signal 
regions with smaller cross sections. 


24 


















would be required to determine the presence of u — t mixing, of which we hnd that 
the channel with two positively charged leptons is an excellent candidate. 

5.3 Same sign top Channel 

Within the framework of SUSY the “smoking gun” signal of the maximally mixed u 
model at the LHC is the final state containing two positive leptons and large missing 
energy. 

The same-sign positive lepton final states are a consequence of the uu —>■ uu 
process, with consecutive decays to txo, or of the ug —)• ug process, where u decay to 
txo and the decay chain of the gluino present one positive top. Though suppressed 
by two powers of the u — t mixing angle, as well as the small branching ratios of 
W to leptons, the -|- Ifij. channel offers a very clean probe of the presence of 
large u — t mixing. Production of u*u* and gu* which would yield two negative sign 
leptons in the final state, contributes only few % to the total signal cross section, due 
to the PDF suppressions. In the context of SUSY, the strong PDF suppression is a 
valuable feature of signals with large ut mixing, as other RP conserving supersym¬ 
metric models can predict same-sign lepton signals with the same amounts of 
and l~l~ events [46]. Furthermore typical RPV models with Baryon violation lead 
to signals with dominant l~l~ [43], while RPV models which also includes lepton 
number violation could lead to dominant [45]. 

Fig. 9 shows the characteristic kinematic distributions of the signal and back¬ 
ground in the -|- Iprp channel. SM backgrounds consist mainly of (P(10) fb level 
rare SM processes (here we consider ttW^tiZ, ZZ,W^W^W~ and W^Z), as well 
as SM ti and W + bb, where the same sign lepton background comes mainly from 
leptonic b decays which yield isolated leptons. With the exception of rare processes, 
the probability that SM processes contain two positive isolated leptons is tiny, yet 
significant due to the large production cross sections. In addition, both the amount 
of missing energy and the transverse mass of the signal events are much larger than 
in the SM backgrounds, as shown in the right panel of Fig. 9. 

In order to maximise the signal significance in the same-sign positive lepton 
channel, here we employ a set of cuts: 

Ni+ {pT > 20 GeV, T] < 2.5) = 2 , > 120 GeV , 

Nj{pT > 30 GeV, p < 2.5) > 1, rriTilt) > , 

W-(pt > 20 GeV, 7] < 2.5) = 0 (5.2) 

where refer to isolated positive/negative leptons and If refers to the highest 
Pt positive lepton in the event. For the purpose of illustration, we set = 

(200,300,400) GeV, for rriu = (450,700,950) GeV respectively. 

® Other non-supersymmetric models such as Z', Composite Top and extra-dimensional models 
can also produce final states with two positive leptons. 


25 




^[GeV] 



rriT I GeV] 



[ GeV] 



rriT [ GeV] 


Figure 9: Kinematic distributions relevant for the signal search in the + Ipj, channel. 
All histograms are normalised to unit area. The label “rare proc.” includes SM production of 
ttW,ttZ, ZZ,W^W'^W~ and WZ. In the top panels, we are showing only the benchmark point 
with (to*, = (450,200) GeV. In the bottom panels, the values in the legend represent different 

{mu,mjgg) benchmark points. 


Our choice of cuts on Ifij, and tut was somewhat optimised to probe high rriu in a 
non-compressed mass spectrum scenario. Again, it is important to note that for the 
purpose of probing the compressed region of rriu, parameter space it is benehcial 
to relax the cuts on missing energy and transverse mass. Bottom panels of Figure 9 
illustrate the point. The benchmark point of = (700,400) GeV displays a 

much softer spectrum of missing energy and transverse mass than the corresponding, 
non-compressed point of = (700,200) GeV, suggesting that a cut on Ifij. 

and rriT lower than the one we suggest in Eq. 5.2 could be more optimal in the 
compressed spectrum scenarios. 

Gontinuing, Table 7 shows an example outflow for two benchmark parameter 
points. The requirement on the lepton multiplicity efficiently reduces the ti and W -|- 
bb backgrounds, while the cut efficiently suppresses the rare-process contribution 
to the total event yield. We hnd that a minimal set of cuts in Eq. 5.2 results in a 
factor of ~ 10 — 15 improvement in S/B, at a cost of ~ 50% in signal efficiency. 

Table 8 shows a summary of results on Eve benchmark points in the rriu^m^g^ 
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rriu = 450 GeV, = 200 GeV, nig = 2 TeV 


-k + jj 

uu 

ug 

tt 

W-kjets 

rare proc. 

S/B 

S/y/B (3000 fb”^) 

Ni+ 

= 2,Ni-=0 

0.21 

0.037 

0.067 

0.022 

19.0 

0.013 

3.1 

Ng 

> 1 

0.20 

0.037 

0.033 

0.022 

14.0 

0.017 

3.4 

IpT 

> 120 GeV 

0.12 

0.034 

0.033 

< 0.01 

1.2 

0.12 

7.6 

niT 

> 200 GeV 

0.080 

0.019 

< 0.01 

< 0.01 

0.45 

0.21 

7.9 



mu 

= 950 

GeV, m 

xo = 400 

GeV, mg = 

2 TeV 


-k + jj 

uu 

ug 

tt 

W-kjets 

rare proc. 

S/B 

o 

o 

o 

1 

Ni+ = 2, Ni- = 0 

0.031 

0.014 

0.067 

0.022 

19.0 

0.0023 

0.56 

Ar,>l 

0.030 

0.014 

0.033 

0.022 

14.0 

0.0031 

0.64 

> 120 GeV 

0.026 

0.014 

0.033 

< 0.01 

1.2 

0.033 

2.0 

rriT > 400 GeV 

0.014 

0.0068 

< 0.01 

< 0.01 

0.081 

0.26 

4.0 


Table 7: Example cutflow in the GG + fj, channel. The entries show cross sections in fb after 
each consecutive cut. The label “rare proc.” includes SM production of tiW,ttZ, ZZ,W^W'^W~ 
and WZ. We compute tot using the hardest positive lepton in the event and fj,. 


m^Jmu ( GeV) 

450 

700 

950 

200 

(0.21, 7.9) 

(0.31, 7.5) 

(0.26, 4.1) 

400 

- 

(0.22 , 5.3) 

(0.26, 4.0) 


Table 8: Summary of reach for benchmark points in the same sign positive lepton channel. The 
table entries show S/B and Sj^/B at 3000 fb”^ respectively which can be achieved at LHC Run 
II in the same-sign positive lepton channel. The masses of squarks and neutralinos are listed on in 
the topmost row and leftmost column respectively. All results assume wig = 2 TeV. 


space. We find that the LHG can achieve a signal signihcance of 5a for u with masses 
of up to ~ 700 GeV for a neutralino mass < 400 GeV. A signihcance higher than 
3a can be achieved for masses up to ~ 1 TeV, assuming the neutralino masses of 
< 400 GeV, suggesting that the same sign positive lepton channel could rule out the 
model up to ~ 1 TeV at the high luminosity LHG. 

The reach of the single top channel at 300 fb“^ is comparable to the reach of the 
same-sign positive leptons search at 3000 fb~^. In case a signal is observed in the 
single top channel at 300 fb~^, the high luminosity LHG should be able to probe and 
measure possible large squark mixings. Gonversely, in case no signal is observed in 
the single top channel at 300 fb~^, the parameter region giving raise to measurable 
same-sign positive leptons yield at High Luminosity LHG will already be ruled out. 
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6 Conclusions 


The canonical paradigm of snpersymmetry typically assumes ultra-violet completions 
with a flavour symmetry, thereby leading to aligned or diagonal soft masses. Yet, 
the possibility of large off diagonal entries represents a viable and interesting option 
to be investigated in the context of the MSSM. 

In this paper we studied an example model of supersymmetry with large mixing 
in the right-handed up-type squark mass matrix in the context of extended gauge 
mediation. Our analysis includes detailed consideration of constraints imposed on 
the parameter space by flavour observables, as well as prospects for LHC Run II 
to discover models with large squark mixings. We hnd that the single-top hnal 
signatures predicted by our model can be accessible at LHC14 with 300 fb“^. A 
more distinctive feature of the MMUT model, the same sign positive lepton hnal 
state which signals a large squark mixing, can be probed at the high luminosity 
LHC. A combination of the two searches could provide useful insight into the havor 
mixing properties of the light squark state. 

As models with large squark mixings typically suffer from an increased degree of 
hne tuning, we take a moment to discuss the naturalness of the MMUT model. The 
parameters determining the tuning in the MSSM are the dimensionful terms entering 
in the corrections of the up type higgs soft mass: 

= 6 Tr {yimQVu + vlmuyu + + ... , 

eventually leading to a hne tuning in the EWSB correction. From this expression, 
one concludes that the relevant entries of the squark mass matrix in the evaluation 
of the tuning are the ones projected onto the Yukawa directions, essentially the (3, 3) 
entries. 

The reference scenario to which we compare the degree of hne tuning is natural 
SUSY [90], where the lightest right handed squark is a pure stop, not mixed in havor. 
The tuning in natural SUSY is set by the LHC bound on right handed stop which, 
in the case of 200 GeV neutralino, is around tuq ^ 650 GeV [73]. In analogy with 
[6], we dehne a parameter measuring the departure from natural SUSY by dividing 
the (3, 3) entry in our MMUT scenario with the minimal one of natural SUSY: 

^ _ Imj + 

— 2 2 2 ’ 

TUq UIq TTIq 

where we labelled with mi and m 2 the lightest and next to lightest up-squark eigen¬ 
state, respectively. 

In the simplihed model of Section 3 we considered only one of the lightest eigen¬ 
states. We found that the lower bound on such state from LHC 8 TeV searches is 
around 700 GeV (neglecting here the possibility of = 450 GeV which is a pecu¬ 
liar very compressed point; see Fig. 7). In order to provide a quantitative estimate 





of we should specify also the value of the other up-squark eigenstate. In Section 
3 we assumed that the next to lightest squark eigenstate is decoupled from LHC 
physics, and taking inspiration from the gauge mediation model of Section 2.2, we 
can consider it to be at the 0(5 — 10) TeV scale. In this case the MMUT model 
would be considerably more tuned than natural SUSY by a factor ^ ~ 30 — 100. In 
the most optimal scenario, instead, we can assume that the same LHC bound on 
the lightest squark eigenstate applies also to the next to lightest eigenstate. In this 
case we apply a common bound of ~ rriu^ > 700 GeV, and the tuning of the 
MMUT model with respect to natural SUSY reduces to ~ 1.2. Hence we hnd that 
in the optimal case, the MMUT model also represents a slightly more un-natural 
SUSY scenario. The reason is that the LHC bound on the lightest squark state in 
the MMUT model is higher than the LHC bound on a pure stop eigenstate, since in 
the former case we have extra production modes (in particular process 1 in Fig. 6; 
see also the cross Section plot of Fig. 7). We conclude that generically the MMUT 
model will be at least slightly more tuned than natural SUSY. 

Future studies of models with large squark mixings would beneht from including 
the possible scenario where the neutralino is not stable on collider scales, adding 
typically two extra displaced photons to the signatures we discussed, which could 
potentially be accessed at the LHC Run H. 

On the model building side, we note that we did not address the issue about the 
dynamical origin of supersymmetry breaking. It would be interesting to explore this 
aspect at greater depth, as well as to evaluate the model’s effective level of tuning 
from a UV perspective. 

Acknowledgments 

We are grateful to Gilad Perez for discussions and insightful comments on a prelim¬ 
inary version of this paper. We also thank David Shih and Jared Evans for useful 
comments on the draft. A.M. would like to thank also Riccardo Argurio, Lorenzo 
Calibbi, Athanasios Dedes, Didar Dobur and Diego Redigolo for interesting discus¬ 
sions. M.B. would like to thank the University of Kansas phenomenology group for 
their hospitality during the hnal stages of the project. M.B. and A.M. are supported 
in part by the Belgian Federal Science Policy Office through the Interuniversity At¬ 
traction Pole P7/37. A. M. is a Pegasus FWO postdoctoral Fellowship. A. M. is also 
supported in part by the Strategic Research Program High Energy Physics and the 
Research Council of the Vrije Universiteit Brussel. 

^°Note that here one should rely on some extra mechanism (such as the NMSSM) to obtain the 
correct Higgs mass. 


29 



A Soft terms from R-symmetric hidden sector 

In this section we obtain the soft term contribution for the model studied in this 
paper, i.e. an hidden sector with a discrete Z4 R-symmetry (2.3) coupled via super¬ 
potential interaction (2.2) to the up type quark. 

The cases of non-R symmetric hidden sectors coupled via messenger matter cou¬ 
pling to the MSSM have been completely classified in [11]. The coupling with the 
same structure we are considering has been denoted as /13 in their classification, 
and it consists of the messenger-matter coupling 

3 

W = Xj2ciUi(jyih (A.l) 

i=l 

where the messenger are assumed to be part of a supersymmetry breaking sector 
coupling to a spurion as 

WsUSYbr = -A(0101 + 0202) (A.2) 


with X = M + e^F. 

In this case the contribution to the soft masses can be extracted from the formulas 
of [11] and results 


^UiFu, = - 




du ,2 A 

16^2^ 


(A.3) 


^ CijX^du ( X'^d^ + X^du - 2 ^ Crgl ) 

r=l,3 


drriQ.. = 


5m|r = 

n.u 


2567r4 

duX^ 
'25(hd 
MuX^ 
' 2567r4 


(yhc.y)ij h? 
Tr(yhc.?/) A^ 


where du = 2 and = 4, Ci = 2/5, (^2 = 0,6*3 = 4, = CiCj and we used that 

4 + 4 + 4 = 1. 

These results have been obtained by studying the threshold corrections to the 
wave function renormalization of the MSSM matter helds induced by integrating out 
the messengers. Precisely, in the case of the up type quark (the other sfermions are 
analogous), the results (A.3) arise from the correction in the Kahler potential 


j rf-9Z„.„,(|X|)C//C/, D + [F(axZ„,)FlU, + A.c.] (A.4) 

leading to 

^UiUj = F{dxZuiUj) (A.5) 

= -\F?{dxdx*ZuiUj) + \F\^{dxZuiu4{dx*Zu,,uj) (A.6) 
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where the second term in the mass squared arises from integrating out the F term 
of the MSSM held. From this procedure it is clear that the contributions to the 
soft masses (at two loop) can be divided into a contribution coming from the second 
derivative of Zu^uj with respect to X and X*, and another one coming from the 
square of the hrst derivative of Zu^uj, be. the A-term squared. Indeed the authors 
of [11] splitted the contributions to schematically into + |Ap to 

make it explicit. In the formulas (A.3) the A-term contribution to the soft masses is 
the second term in the big round parenthesis. 

The last term in in (A.3) is coming from one-loop correction to the soft 

mass, it is not at leading order m F/M but plays a crucial role in lowering the squark 
mass. The precise expression for the h function is 

h{x) = — 2) ln(l — x) — {x + 2) ln(l -|- x) j (A.7) 

and it is an even function of x. 

Now, in order to obtain the complete set of soft terms induced in our model 
by the R-symmetric SUSY breaking sector (see eqn. (2.3)), we adopt the following 
strategy. Consider a double copy of the above non R-symmetric model, with two 
different spurious Xi and X 2 , but with the same coupling to the MSSM matter held 

3 

W = X CiUi{(l)i(j)2 + 0304) + Yi(0i0i -(- 0202) + ^2(0303 + 0404) (A.8) 

i=l 

As in the case above, the one loop corrections to the A-terms and the two loop 
corrections to the soft masses are encoded into the wave function renormalization for 
the MSSM held which receive in this case two additive contributions 


ZuiUj — ZuiuA\Xi\) + (IA2I) 


(A.9) 


Note that the functional form of Z^.^Al^il) cind Zu^Ujil^A) is the same. 
The correction induced by the Kahler potential are then 


d‘0Z^„UlU, D + FA,Z„,^,(\X 2 \))fIUj + h.c. 


+ (\Fi\\dxA-,Z„,(\X,\) + \F2\‘(d>c,dx-,Z„,(\X^\)\UlUi 


Now, in the special case in which 


Xi = M + e^F X 2 = M-e‘^F (A.IO) 

^^The two sectors are coupled only via the up type quarks and the correction of one sector to 
the other are loop suppressed and negligible, unless there is a huge hierarchy between Xi and X 2 . 
This is analogous to the situation for two SUSY breaking hidden sector in gauge mediation [91]. 
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we have that the terms with hrst derivatives of cancel out and we are left only 
with the second derivatives contribution to the soft masses, which are two copies of 
the same expression. Hence for this model, with the particular choice (A. 10), the 
induced soft terms are twice the following contributions 

I - 2 Crgl J (A. 11) 

\ r=l,3 / 

In the special choice (A.10), we can rotate the messenger helds in a new basis such 
that the theory (A.8) is equivalent to the following model 

3 

A ^ ^ CjI7j(0a0bA0c0(i)“l“Af (0(j0Q,-|-0b0i)A0c0c“l“0(i0d)“l“H(0Q0j-(-0j0Q,-|-0(;0(i'h0(i0c) 

(A.12) 

where we denoted Y = 9‘^F, which is exactly two copies of the R-symmetric model 
considered in the paper, see eq (2.3) and (2.2). We hence conclude that the soft 
terms induced by the R-symmetric model at leading order in F/M are given by 
(A. 11). Moreover we also computed explicitly the one loop corrections (suppressed 
in F/M) to the soft masses in the R-symmetric model, and found that they are 
equivalent to the one induced by the non R-symmetric model, i.e. the last term in 

in (A-3). 

Hence the total contribution to the soft terms for the model discussed in the 
main body of the paper, i.e. (2.3) and (2.2), is indeed the one quoted in (2.5). We 
stress that no A-terms are generated, and the |Ap term is not present in the two 
loop corrections to the soft masses. 
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